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Abstract. We compare the fractions of the hard and geometric cross sections 
as a function of impact parameter. For a given definition of central collisions, 
we calculate the corresponding impact parameter and the fraction of the hard 
cross section contained within this cut. We use charm quark production as a 
definite example. 

In this note, revised from Ref. standard nuclear density distributions are 
described and the resulting geometrical overlap in nuclear collisions is calculated. 
We then compare the hard process cross section to the total geometric cross section 
as a function of impact parameter and discuss how they are related to the collision 
centrality. 

A three-parameter Woods-Saxon shape is used to describe the nuclear density 
distribution, 

l+Lo(r/R A ) 2 

PA(r) = p — p~~r7~r ■> I 1 

1 + exp((r - Ra)/z) 

where Ra is the nuclear radius, z is the surface thickness, and u> allows for central 
irregularities. The electron scattering data of Ref. Q is used where available for 
Ra-, z, and lo. When data is unavailable, the parameters u> = 0, z = 0.54 fm 
and R A = 1.19A 1 / 3 - 1.61A- 1 / 3 fm are used. The central density po is found 
from the normalization J d 3 rpA(r) — A. For results with other nuclear shape 
parameterizations, see the appendix of Ref. ||. The nuclear shape parameters 
are given m Table 0. See also Ref. {§ for more detailed discussion of nuclear density 
distributions. 

In minimum bias (impact parameter averaged) AB collisions we expect the 
production cross section for hard processes to grow approximately as 

hard hard / a r>\ct /n\ 

CT AB =°pp \ AB ) > ( 2 ) 
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A 


Ra (fm) 


z (fm) 


u> 


Po (fm" 3 ) 


16 


2.608 


0.513 


-0.051 


0.1654 


27 


3.07 


0.519 


0. 


0.1739 


40 


3.766 


0.586 


-0.161 


0.1699 


63 


4.214 


0.586 


0. 


0.1701 


110 


5.33 


0.535 


0. 


0.1577 


197 


6.38 


0.535 


0. 


0.1693 


208 


6.624 


0.549 


0. 


0.1600 



Table 1. Nuclear shape parameters taken from Ref. [0. 



where a = 1 when no nuclear effects are included. However, central collisions are of 
the greatest interest since it is there that high energy density effects are most likely 
to appear. Central collisions contribute larger than average values of Et to the 
system, in the 'tail' of the Et distribution, da/dEx- We would like to determine 
which impact parameters are important in the high Et tail, i.e. what range of b 
may be considered central. We now define the central fraction of the hard cross 
section, Eq. (Q), and the central fraction of the geometric cross section. We then 
discuss how the two are related. 

Considering only geometry with no nuclear effects, a = 1 in Eq. (|[), the inclu- 
sive production cross section of hard probes increases as 

^AB d = <J^ d T AB (b)d 2 b (3) 

and the average number of hard probes produced at impact parameter b is N AB (b) = 
<7pp ld T AB (b) where T AB (b) is the nuclear overlap integral, 

Tab (b) = J d 2 sT A (s) T B (b - s) , (4) 

and Ta = / dzp A (z, s) is the nuclear profile function. The nuclear overlap functions 
for Pb+Pb and Au+Au collisions are shown in Fig. [j] as a function of impact 
parameter. Integrating Tab over all impact parameters we find 

J d 2 bT AB (b) = AB . (5) 

The central fraction /ab, equivalent to the fraction of the total hard cross 
section, is defined as 



/ab = ^ J bdbT AB (b) , 



(0) 
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0.0 2.5 5.0 7.5 10.0 12.5 15.0 

b (fm) 

Fig. 1. The nuclear overlap function Tab(^) as a function of impact parameter b 
for Pb+Pb and Au+Au collisions. 

where b c is the central impact parameter and b < b c are central. To make a similar 
'centrality cut' in pA collisions, the fraction 

/a = ^ J'bdbT A (b) (7) 

would be used. Fig. |], taken from Ref. shows the increase of Jab with b c for 
several symmetric, AA, systems. Note that /aa ~ 1 when b c sa 2Ra- For example, 
if we choose cr cen trai = 0.1(TAB d ; this corresponds to b c = 2.05 fm in Au+Au collisions 
and b c = 1.05 fm in O+O collisions. If we instead chose cr ccntra i = 0.01<7 A g d then 
b c = 0.52 fm in Au+Au and b c = 0.33 fm in O+O collisions. 

Note however that /ab is not the fraction of the geometric cross section which 
includes both hard and soft contributions. The geometric cross section in central 
collisions is found by integrating the interaction probability over impact parameter 
up to b c , 

Vgeo(b c ) = 27r / bdb [1 - exp(-TABCjvAr)] ■ ( 8 ) 
Jo 

The nucleon-nucleon inelastic cross section, <tnn, is ~ 32 mb at SPS energies and 
grows with energy. It is expected to be ~ 60 mb at LHC energies. The fraction of 
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the geometric cross section is 



Jgeo — 



'gco 



(be) 



(9) 



'gco 



In central collisions, where Tab is large, the impact parameter dependence is simple, 
Ogco(frc) oc b 2 c . However, in peripheral collisions where the nuclear overlap becomes 
small, cr geo (b c ) deviates from the trivial b 2 . scaling. Deviations from this scaling do 
not occur until b c w 2Ra in symmetric systems. 
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Fig. 2. The central fraction of the hard cross section as a function of impact 
parameter cut b c for several symmetric systems. 

Figure || shows the numerical result, Eq. (||), relative to the integral where 
b c — > oo, for the same systems as in Fig. We have used &nn = 32 mb in Eq. (|^). 
A negligible difference in the most peripheral collisions can be expected if 60 mb 
were used instead. The growth of the fraction of the geometric cross section is 
slower than that of the hard fraction, /ab- Indeed at b c « 2Ra, / gco ~ 0.75. 

The total geometrical cross section for a variety of colliding nuclei is given 
in Table 0. We have also calculated the impact parameter b c corresponding to 
/gco = 0.05, 0.1, and 0.2 or the central 5%, 10% and 20% of all collisions. The 
impact parameter corresponding to / gco = 0.2 is b c ss l.Q4R^ when symmetric 
systems are considered. In asymmetric collisions, b c < Ra when / gco = 0.2. If 
smaller centrality cuts are imposed, the impact parameters are reduced by factors 
of \/2 and 2 for / geo = 0.1 and 0.05 respectively. Note that when B ^> A, in the 
smaller nucleus is embedded in the larger with a 10% or smaller centrality cut. 
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b c (frn) 




A + B 


(Tgco (b) 


/ g eo = 0.05 


/ g eo = 0.1 


/goo 0.2 


16+16 


1.18 


1.37 


1.94 


2.74 


16+27 


1.45 


1.52 


2.15 


3.04 


16+40 


1.75 


1.67 


2.36 


3.34 


16+63 


2.15 


1.85 


2.62 


3.70 


16+110 


2.69 


2.07 


2.93 


4.14 


16+197 


3.42 


2.33 


3.30 


4.66 


16+208 


3.59 


2.39 


3.38 


4.78 


27+27 


1.76 


1.67 


2.37 


3.35 


27+40 


2.09 


1.82 


2.58 


3.65 


27+63 


2.53 


2.01 


2.84 


4.01 


27+110 


O 1 1 

3.11 


2.22 


3.15 


A A C 

4.45 


27+197 


3.89 


2.49 


3.52 


4.98 


27+208 


4.08 


2.55 


3.60 


5.09 


40+40 


2.45 


1.98 


2.79 


3.95 


40+63 


2.93 


2.16 


3.05 


4.32 


40+110 


3.55 


2.38 


3.36 


4.76 


40+197 


4.38 


2.64 


3.73 


5.28 


40+208 


4.58 


2.70 


3.82 


5.40 


63+63 


3.46 


2.34 


3.32 


4.69 


63+110 


4.14 


2.56 


3.63 


5.13 


63+197 


5.04 


2.83 


4.00 


5.66 


63+208 


5.25 


2.89 


4.09 


5.78 


110+110 


4.86 


2.78 


3.93 


5.56 


110+197 


5.82 


3.04 


4.30 


6.09 


110+208 


6.06 


3.10 


4.39 


6.21 


197+197 


6.88 


3.31 


4.68 


6.62 


197+208 


7.13 


3.37 


4.76 


6.74 


208+208 


7.39 


3.43 


4.85 


6.86 



Table 2. Values of the geometric cross section and the impact parameter at which 
/goo = 0.05, 0.1 and 0.2 respectively for several colliding systems. 
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b c (fm) 

Fig. 3. The fraction of the total geometrical cross section as a function of impact 
parameter cut b c for several symmetric systems. From left to right the curves are 
16+16, 27+27, 63+63, 110+110, and 197+197. 



For the same systems, Table g gives the value of the nuclear overlap at b = 0, 
Iab(O), and the fraction of the hard cross section, /ab, corresponding to 5%, 10%, 
and 20% of the geometric cross section, calculated with the b c values given in 
Table ||. For example, the central 10% of the total geometric cross section is obtained 
when b c ~ 4.7 fm in Au+Au collisions, more than twice the corresponding impact 
parameter for the same percentage of the hard cross section. In this case, 10% of 
the geometric cross section encompasses « 40% of the hard cross section. In fact, 
a 10% cut on the geometric cross section corresponds to 30-40% of the hard cross 
section for all systems considered. Even the central 5% of collisions encompasses 
17-23% of the hard cross section while a less stringent cut of 20% garners 52-66% 
of all hard probes produced before nuclear effects are considered. 

In Fig. [|, we show the ratio of /ab relative to the geometric ratio for the same 
systems as in Figs. |2| and 0. The hard fraction grows more slowly relative to the 
geometric fraction in smaller systems, 16+16 and 27+27, but otherwise the results 
are similar. Figure || shows the relative ratios for the asymmetric systems 16+197, 
27+197, 63+197, 110+197 and 197+197. In this case, the relative ratios cluster 
even closer together than for symmetric systems. Thus the larger nucleus sets the 
scale for both the hard and geometric cross sections in asymmetric systems. 

The most appropriate way to obtain the number of hard probes produced in a 
central collision is to calculate b c from the geometric cross section and then, with 
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/ab 




A + B 


Tab(0) (mb r ) 


/ g eo = 0.05 


/geo = 0.1 


/goo 0.2 


16+16 


0.79 


0.166 


0.306 


0.522 


16+27 


1.16 


0.178 


0.326 


0.550 


16+40 


1.44 


0.180 


0.330 


0.558 


16+63 


1.89 


0.185 


0.340 


0.576 


16+110 


2.41 


0.174 


0.325 


0.565 


16+197 


3.22 


0.167 


0.316 


0.560 


16+208 


3.18 


0.164 


0.312 


0.555 


27+27 


1.75 


0.191 


0.348 


0.580 


27+40 


2.23 


0.195 


0.355 


0.591 


27+63 


3.00 


0.201 


0.365 


0.606 


27+110 


3.94 


0.191 


A O C A 

0.354 


0.599 


27+197 


5.32 


0.184 


0.344 


0.595 


27+208 


5.27 


0.181 


0.340 


0.591 


40+40 


2.92 


0.202 


0.365 


0.604 


40+63 


4.05 


0.210 


0.379 


0.623 


40+110 


5.53 


0.204 


0.371 


0.619 


40+197 


7.64 


0.198 


0.366 


0.618 


40+208 


7.57 


0.195 


0.361 


0.613 


63+63 


5.77 


0.220 


0.395 


0.642 


63+110 


8.19 


0.218 


0.393 


0.643 


63+197 


11.61 


0.213 


0.389 


0.644 


63+208 


11.57 


0.211 


0.385 


0.640 


110+110 


12.43 


0.222 


0.398 


0.648 


110+197 


18.54 


0.223 


0.400 


0.652 


110+208 


18.64 


0.221 


0.397 


0.649 


197+197 


29.32 


0.229 


0.410 


0.663 


197+208 


29.82 


0.229 


0.409 


0.662 


208+208 


30.42 


0.229 


0.409 


0.664 



Table 3. Values of Tab (0) and the fraction of the hard cross section for / gco = 0.05, 
0.1 and 0.2 respectively in several colliding systems. 
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geo 



Fig. 4. The fraction of the hard cross section as a function of the total geometrical 
cross section for the symmetric systems shown in Figs. 2 and 3. From left to right, 
the curves are 197+197, 110+110, 63+63, 27+27, and 16+16. 

this & c , calculate /ab with Eq. (|^). Assuming no nuclear effects on the hard probe 
production cross section, the average hard process rate at b = in Au+Au collisions 
is 

-hard . 



N 



AL! 



(0) 



hard 



Tab(O) 



(10) 



where cr^ d is the total hard process production cross section in pp interactions. 
The rate in the impact parameter interval < b < b c is the ratio of the hard to 
geometric cross sections integrated over b, 



using Eqs. 



-^AB (be) = 

|), (§, and (| 

R h AB d (bc) 



[ b <= J^-hard hard a r> f 

Jn aa AB a pp At>L 



/AB 



^geo(^c) CTgeo fgco 

In Fig. ||, the ratio 
7V A a B rd (6 c ) 1 ABf AB 



-hard 



(ii) 



(12) 



is shown for the same set of symmetric systems as in Figs. 2-4 as a function of b c . 

As a specific example, the average number of cc pairs produced at b = in 
Au+Au collisions is 



Ar AB (o) = <C T A B (o) 



(13) 
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At */s = 200 GeV, with MRS D-' parton distributions, a™ = 0.344 mb @ and 
7ab(0) = 29.3/mb, resulting in « 10 cc pairs per Au+Au collision at b = 0. With a 
10% centrality cut, the average number of cc pairs produced in the range < b < b c 
is 

Nf M = ^^. (14) 

Since the central 10% of the geometric cross section corresponds to 40% of the hard 
cross section, there are ss 7.7 cc pairs in the 10% most central events. 
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